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The Planck satellite has identified several patches of sky with low polarized dust emission, obvious 
targets for searches for the cosmic-microwave-background (CMB) B-mode signal from inflationary 
gravitational waves. Still, given the Planck measurement uncertainties, the polarized dust fore¬ 
grounds in these different candidate patches may differ by an order of magnitude or more. Here 
we show that a brief initial experiment to map these candidate patches more deeply at a single 
high frequency can efficiently zero in on the cleanest patch(es) and thus improve significantly the 
sensitivity of subsequent B-mode searches. A ground-based experiment with current detector tech¬ 
nology operating at > 220 GHz for 3 months can efficiently identify a low-dust-amplitude patch 
and thus improve by up to a factor 2 or 3 on the sensitivity to cosmic B modes of the subsequent 
lower-frequency deep integration. A balloon experiment with current detector sensitivities covering 
the set of patches and operating at ~ 350 GHz can reach a similar result in less than two weeks. This 
strategy may prove crucial in accessing the smallest gravitational-wave signals possible in large-held 
inhation. The high-frequency data from this exploratory experiment should also provide valuable 
foreground templates to subsequent experiments that integrate on any of the candidate patches 
explored. 

PACS numbers: 


A significant effort is now underway to detect the 
B-mode signal [1] induced by inflationary gravitational 
waves [2] in the polarization of the cosmic microwave 
background (CMB). Any such search must, however, dis¬ 
tinguish the cosmic signal from the contamination by 
Galactic foregrounds, in particular polarized emission 
from dust, as recent experience has shown [3H6]. Recent 
data from the Planck satellite on polarized dust emis¬ 
sion at 353 GHz has shown that dust polarization may 
be higher than had been expected, even in low-intensity 
regions at high Galactic latitudes [7]. The data also indi¬ 
cate that the amplitude of dust-intensity fluctuations in 
a given patch of sky is not necessarily a reliable proxy for 
the dust-polarization fluctuations, as had been assumed 
before. 

The importance of the detection of inflationary B 
modes for cosmology motivates us to reconsider, in the 
light of this new information, every possible avenue to 
reduce the contamination from dust in GMB B-mode 
maps and thus maximize the sensitivity of B-mode ex¬ 
periments. An unambiguous detection will of course re¬ 
quire multifrequency component separation, and a pu¬ 
tative cosmic signal must have the correct statistical 
properties [5]. Still, several experiments that focus on 
a small 400 square degrees) patch of sky [9] (e.g., BI- 

GEP3, PolarBear.), to target the recombination peak 

in the B-mode power spectrum, will do better if they 
integrate on cleaner (i.e., lower dust-polarization fore¬ 
ground) patches of sky. An obvious question for these 
experiments is which of the ^ 100 ^400-square-degree 
patches of the sky to survey. In the past, this selection 
was based on available information on dust-intensity fluc¬ 
tuations. Now that we have 353-GHz polarization maps 


from Planck [3, the selection will be based largely upon 
that information. 

The silver lining for B-mode searches of the Planck 
dust-polarization maps is the identification of a handful 
of patches that may potentially have much lower dust- 
polarization amplitudes than the average patch. The cru¬ 
cial point, though, is that due to the limited Planck sen¬ 
sitivity, the measurements in these low-dust-polarization 
patches are noise-dominated, and it is impossible to de¬ 
termine at this point which of these candidate patches is 
cleanest of dust-polarization foregrounds. The dust am¬ 
plitudes could plausibly vary by an order of magnitude 
or even more. 

Small-sky experiments are hence faced with a difficult 
choice: which of Planck’s best patches should be chosen 
for the deep integration required to optimize the sensitiv¬ 
ity to B modes? Should they choose one at random and 
hope for the best? or perhaps split the total observation 
between multiple patches? 

An alternative strategy, which we consider here, is to 
conduct a brief initial exploratory survey, at a high (dust- 
dominated) frequency, of the few relatively clean Planck 
patches to identify the cleanest of them. As we will show, 
a B-mode experiment that then focuses on this cleanest 
patch may potentially be ~ 20% — 66% more sensitive 
than it would be without this initial stage of exploration. 
More precisely, we use as a figure of merit the lowest up¬ 
per bound achievable by a null experiment on the tensor- 
to-scalar ratio r. We then compare the bound obtained 
in several different scenarios for the exploration phase 
to that obtained with no exploration. We show that the 
~ 20%—66% improvement can be achieved with an initial 
3-month ground-based survey at 220 GHz, with current 



2 


detector technology, or with a two-week balloon experi¬ 
ment at 353 GHz. We also consider briefly several strate¬ 
gies to optimize the initial stage of exploration. We will 
see that the proposed initial exploration may allow exper¬ 
iments to more efficiently access the lowest gravitational- 
wave amplitudes expected in large-held models [iniiii]. 

We use standard techniques (see, e.g.. Ref. [H]) to esti¬ 
mate the sensitivity of a given experiment to B modes in 
the presence of foregrounds and lensing. The Fisher fore¬ 
cast for the error in the measurement of the amplitude A 
of a power spectrum Ce, with some assumed multipole- 
moment ^ dependence, is mm. 



We assume that the likelihood function is Gaussian in the 
vicinity of its maximum To emphasize the rough 

nature of this approximation we use quotation marks 
around “la” when referring to the error aA in the am¬ 
plitude A. 

We now consider an experiment to identify through 
measurements at a single high (dust-dominated) fre¬ 
quency r'dust the cleanest of several candidate patches 
of sky. We assume that the observation time is divided 
into several steps, each dedicated to a different patch of 
sky. We then need to estimate the dust amplitude A in 
a targeted patch within the observation time allotted to 
a single step of the exploration stage. For a given sky 
coverage /sky, the “Itr” error to the value of A estimated 
from an individual multipole moment £ is [anil mi, 



where A denotes the estimated value of A, is the 
lensing B-mode contribution and w“^(tstep) is the inverse 
weight per solid angle given an observation time tstep- 
The quantity 1 — a parametrizes the level of de-lensing 
[IHHini that is applied to the data. To be conservative, 
we assume a = 1 (no de-lensing) unless stated otherwise. 
The total “Icr” error in the measurement of A over a time 
Gtep is thus. 


a 


A 


'/sky 'y {2£ + l){Cf)^ 

. ^ tliCt + 
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where Cf = Cf /A = 27rt'“'"/[£(f-|- 1)] encodes the I de¬ 
pendence of the dust power spectrum (which we assume 
to be the same across the sky, as suggested by Ref. 0 ) 
and £min = 40 is the lowest multipole included in the 
analysis. We emphasize that we are interested in the dust 
amplitude A in the particular patch we are considering 
for prolonged observation, and hence no cosmic-variance 
term appears in Eq. ([^. 


After the optimal patch is identified through explo¬ 
ration, a prolonged stage of integration is performed over 
the chosen patch at a GMB-domimated frequency vcmb 
(and as we discuss below, perhaps additional frequen¬ 
cies). To estimate the gain achieved by exploration, we 
compare the sensitivity to primordial B modes from the 
integration. The smallest amplitude in a patch p de¬ 
tectable at “Icr” (with a total observation time T spent 
on the chosen patch), according to Eq. Q, is then [T^ 


/sky 

( yJ{2l+\)Cf \ 

2 ^ 

i^QC/-|-/3Ap(l/CMB)C'/+/sky«.(T)-le''‘''b ) 


(4) 

where Ap(i/CMB) is the dust amplitude in the patch p 
(extrapolated to frequency jccmb), Cf denotes the £ de¬ 
pendence of the IGW B-mode power spectrum and j3 is 
the fraction of the foreground left in the vcmb map, after 
component separation is performed (in a multifrequency 
experiment); /? = 1 with only one vcmb- We set the 
sample variance of the primordial signal to zero in this 
expression (to compare with the null hypothesis). 

An important issue if that of the bias stemming from 
dust power in cases where it is not cleaned (e.g. in a 
single-frequency experiment). One could attempt to re¬ 
solve this in a number of ways. First, as discussed in 
|23) . it is possible to fit simultaneously for r and A on 
the integration data, using the different ^-dependence of 
dust and gravitational waves. Another way is to em¬ 
ploy the exploration measurements on the same patch of 
sky, either for standard component separation, or as a 
means of inferring the dust amplitude which could then 
be extrapolated to the integration frequencies and used 
to correct the data. Naturally, full component separation 
with multi-frequency data is the superior method, and we 
will assume that is possible for most of the experimental 
setups we consider below. 

The instrumental noise in a CMB-polarization exper¬ 
iment is determined by the detector-array sensitivity 
s = Sdet /VNdet (where Sdet is the noise-equivalent tem¬ 
perature NET of each detector and A^det is the num¬ 
ber of detectors in the array), the angular resolution 
the sky coverage /sky, and the total observation 
time T (which is reduced in practice by the observing 
efficiency). The pixel noise CTpix is then determined by 
CTpix = s/V^pix, where tpix = T/Npi^ is the observation 
time dedicated to each pixel. Defining the inverse weight 
w~^{T) = Airs'^jT per solid angle, the angular power 
spectrum of the instrumental noise, assuming the experi¬ 
mental beam is approximately Gaussian in shape, is given 
by [21] 


CN ^ ^ 


Npix 




e" = fs^yW-\T)e 




( 5 ) 


where H = 47r/sky and al = 6>^^ij^/(8In2). 
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In Table we list the parameters of the instruments 
we consider for the exploration experiment and the sub¬ 
sequent B-mode integration, where for all experiments we 
fix the NET per detector to be Sdet = 480 /rK-^s (which 
is already surpassed by most current generation instru¬ 
ments). We note that these parameters are similar to 
those considered, e.g., in Ref. [ 12 ]. Since dust is bright 
at high frequencies, current-generation ground-based or 
balloon experiments with N^et = O(IO^) will most likely 
suffice. For the prolonged integrations we also consider 
more sensitive experiments (see Ref. [H]), with iVdet as 
high as 10®. We set the minimum angular resolution to 
^fwhm = 30' for all experiments, so that the primordial 
peak at £ ~ 80 may be resolved.Optimistic de-lensing 
prospects such as we consider below will require higher 
resolution E-mode polarization data along with decent 
lensing potential estimates (e.g. via tracers such as the 
CIB) [131 [T4|. Finally, we assume an observing efficiency 
of 20% and a sky coverage of /sky = 1%- 



Type 

iVdet 

^obs [GHz] 

iobs [days] 

Explorationl 

Ground 

10^ 

220 

90 

Exploration2 

Balloon 

500 

353 

12 

Integration! 

Ground 

10^ 

150 

1000 

Integration2 

Ground 

10-* 

{90, 150, 220} 

1000 

Integrations 

Ground 

10® 

(90, 150, 220} 

1000 


TABLE I: Table of experimental parameters for the explo¬ 
ration and integration experiments under consideration. 


The Planck experiment reported in [7] the results of an 
analysis of its 353 GHz maps at high Galactic latitudes 
(|6| > 35°), where 352 patches of 400deg^ were analyzed 
in the multipole range 40 < £ < 370, to determine the 
amplitude of dust polarization. This analysis identified 
a group of candidate patches with low dust amplitudes, 
although these amplitudes had large uncertainties. 

In Table |n| we list the amplitudes of Planck’s six best 
candidates, normalized to the primordial gravitational- 
wave B-mode power spectrum at £ = 80 and extrapolated 
down from 353 GHz to 150 GHz, assuming a modihed 
blackbody spectrum with spectral index Pd = 1.59 and 
Td = 19.6Rr, applying the unit conversion factors and 
color correction coefficients as explained in [22]. 


Given the large uncertainties in the measurements of 
the dust amplitudes in these six patches, we cannot claim 
to know how much better one can do in terms of B-mode 
detection by locating the cleanest patch of the sky. We 
can, however, study the potential improvements obtained 
from locating the cleanest patch of sky by simulating 
skies consistent with Planck measurements. 

To do so, we generate an ensemble of simulated skies 
where the dust amplitudes in the candidate set are drawn 
from normal distributions with the best-fit mean and 
variance reported by Planck N(rd,<Jri) (for each patch 
in each simulation, we repeat the draw until a positive 
amplitude is drawn), and calculate ensemble averages for 
all desired quantities. We then use the conversion con¬ 
vention described above to transform between 150 GHz 
and either 220 GHz or 353 GHz. 

In Table HI we calculate the lowest detectable tensor- 
to-scalar ratio for the worst, mean, and best foreground 
amplitude, averaged over the ensemble of simulations, as¬ 
suming a single-frequency experiment (where /3 = 1) as 
well as multifrequency experiments, where after compo¬ 
nent separation a residual of either P = 10% or /3 = 1% 
is left, with and without de-lensing (we consider either 
no de-lensing, a = 1, or an optimistic 90% de-lensing, 
a = 0.1). The parameters for these B-mode-integration 
experiments are taken from Table |T] 

We should mention that there is some redundancy in 
the scenarios covered in Tables HI and IHII With no de- 
lensing, for example, the third integration experiment 
will be lensing-noise limited in a small patch of sky and 
will not exceed the sensitivity of the second experiment, 
unless a larger area of sky is covered, to reduce the 
sample-variance error. For simplicity, we neglect such op¬ 
timizations here as our main concern is the effect of fore¬ 
ground reduction (by means of avoidance) and our hgure- 
of-merit is the lowest bound achievable on the tensor-to- 
scalar ratio under the null hypothesis. Naturally, these 
considerations should play a dominant role when plan¬ 
ning an actual experimental strategy for a given instru¬ 
ment (along with various other considerations, e.g. con¬ 
straints on the observing efficiency in different regions of 
the sky). 

Glearly, there is room for improvement if cleaner 
patches can be identified before the prolonged integra¬ 
tion is done. We therefore turn to investigate how best 
to identify the cleanest patch. 


Td 0.053 0.027 -0.062 -0.020 0.057 -0.031 

ad 0.096 0.098 0.052 0.127 0.122 0.121 


TABLE II: The best-fit dust amplitude Vd and its uncertainty 
Od for the six best patches identified by Planck at 353 GHz. 
The values shown here are normalized to the amplitude of 
primordial B-modes at £ = 80 and extrapolated to 150 GHz 
assuming a modified blackbody spectrum |22 |. 


Given a set of candidate patches for which we have 
rough dust-amplitude estimates, we can choose to forgo 
exploration (and just choose the best candidate accord¬ 
ing to the available data, uncertain as it may be) or to 
pursue a uniform exploration of all patches so as to lower 
the measurement uncertainty in the dust amplitude and 
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d 

Worst Mean 

Best 

Worst 
a = 0.1 

Mean 

0 = 0.1 

Best 

Q = 0.1 

Integration! 

1 

0.034 

0.022 

0.011 

0.03 

0.018 

0.009 

Integration2 

1 

0.025 

0.014 

0.006 

0.019 

0.01 

0.003 

Integration2 

10% 

0.007 

0.005 

0.004 

0.0033 

0.0022 

0.0011 

Integrations 

1 

0.023 

0.013 

0.005 

0.016 

0.008 

0.002 

Integrations 

10% 

0.006 

0.005 

0.003 

0.0024 

0.0014 

0.0006 

Integrations 

1% 

0.0033 0.00310.0029 

0.0007 

0.0005 

0.0004 


TABLE III: Table of expected lower bounds on r with differ¬ 
ent integrations, with no exploration. 


thus allow a more substantiated choice. Alternatively, 
one could adopt an adaptive strategy, whereby the al¬ 
location of exploration time is determined on-the-fly ac¬ 
cording to the gathered results. In this case, patches 
that acquired data indicate are most likely inferior are 
discarded early on, leaving more time to explore more 
deeply the remaining patches and thus zero in efficiently 
on the best ones. 

In Ref. [ 53 , we investigated a variety of adaptive- 
survey strategies for ground-based single-frequency ex¬ 
periments (similar to BICEP2). We showed how to 
construct adaptive survey strategies, based on machine¬ 
learning algorithms, to maximize the sensitivity to pri¬ 
mordial B-modes by spending more time observing lower- 
foreground patches. We concluded that a wisely chosen 
adaptive-survey strategy could improve the upper limit 
to B modes (assuming a null experiment) by a factor of 
2-3. 

Here, though, we consider a limited fixed-time explo¬ 
ration at a dust-dominated high frequency which is then 
followed by a prolonged integration at CMB frequencies. 
In this case, the advantage of adaptive strategies is lim¬ 
ited, first because the measurements are easier and sec¬ 
ondly because there is no price paid in terms of inte- 
gation time (the integration is done separately, by an 
independent CMB-frequency measurement). Neverthe¬ 
less, we compare the performance of the adaptive Up¬ 
per Confidence Bound (UCB) adaptive-survey algorithm 
(which was found in Ref. [53] to be quite effective) to 
a more naive uniform exploration to study the improve¬ 
ments that can be made if a given experiment can move 
quickly between patches at low cost. 

In Fig.[^ we track the average Icr bound on the tensor- 
to-scalar ratio r that may be reached by a subsequent 
integration experiment (we plot the results for the in¬ 
strument Integration2 with (3 = 10% and a = 0.1), given 
each additional day of exploration. At each point in time, 
we plot the achievable bound on the tensor-to-scalar ra¬ 
tio r assuming it observes the patch with the lowest- 
measured dust amplitude by the exploration experiment 
up to that stage. In the top panel we show results from 
a ground-based experiment and in the bottom from a 
balloon instrument, with parameters from Table [ij As 


explained above, we compare three strategies. The black 
solid line shows the average daily bound assuming a patch 
is chosen at random and no exploration performed (hence 
no improvement). Its range of performance is shown in 
grey. The red solid line shows the daily bound as the 
patches are consecutively observed, each for a period of 
15 days in the case of the ground-based instrument (or 
two days for the balloon). During the first observation 
period, there would be no gain compared to choosing at 
random, as expected. Then, as each additional patch 
is observed, the projected bound on r improves as more 
information is gathered to allow an intelligent choice of 
patch for integration. The range of performance of this 
method is shown in red. Finally, the blue solid line shows 
the daily bound with the UCB method, which converges 
onto cleaner patches much more quickly, but requires a 
flexible experiment which can switch between patches ev¬ 
ery three days (or a single day, for the balloon experi¬ 
ment). As can be seen, the worst case when adopting 
this strategy is almost always better than with the other 
methods. 

How well will the B-mode searches (the integrations) 
perform if they are preceded by a stage of exploration 
as described above? In Table |IV| we present the average 
and worst-case B-mode limits achievable by the same in¬ 
tegration experiments considered in Table |HI[ after the 
ground-based exploration experiment in Table has been 
performed to identify the optimal patch for observation in 
each simulation. (The balloon-borne exploration yields 
very similar results, so we do not show them separately). 
We see that the improvement compared to the mean 
without exploration ranges from 20% to 66% in almost all 
experiments, depending on the level of foreground sub¬ 
traction and the amount of de-lensing. Another advan¬ 
tage of the exploration stage is that the worst case is now 
greatly reduced compared to the non-exploration case, as 
can be seen clearly. 

The only case in which there is only minor improve¬ 
ment, < 5%, is with A^det = 10®, 1% foreground resid¬ 
uals, and no de-lensing. There is no improvement here 
as both the foreground residuals and instrumental noise 
are lower than the lensing B modes for all six patches 
and hence there is no gain by focusing on cleaner ones. 
Such a tactic, though, cannot reach the value r ~ 0.002 
that is targeted to fully probe the parameter space for 
large-field-inflation. 


To conclude, we have considered the benefits for B- 
mode searches of an initial exploratory experiment to 
measure at a single high (dust-dominated) frequency the 
dust-polarization amplitudes in several candidate low- 
dust-polarization patches identified by Planck. We have 
shown that the sensitivity of the subsequent B-mode 
searches can be improved, with the identification of the 
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FIG. 1: The projected “Icr” bound on the tensor-to-scalar 
ratio r with the instrument Integration2 with /3 = 10% and 
a = 0.1 (see Table |^, achievable with the information ob¬ 
tained from X days of exploration. We plot the results given 
a ground-based experiment {Top), or a balloon-borne one 
{Bottom), both with parameters taken from Table Three 
methods are compared: no exploration (solid black line); con¬ 
secutive equal-time observation (solid red line); an adaptive 
method (UCB) (solid blue line). The ranges of performance 
are also shown (in corresponding color shades) to demonstrate 
the worst case for each method. After all patches are observed 
for l/6th of the total exploration time, the uniform method 
reaches the same bound as the adaptive one. 



p 

Worst 

Mean (A(%)) 

Worst 
a = 0.1 

Mean A(%)) 
« = 0.1 

Integration! 

1 

0.025 

0.012 (39%) 

0.0213 

0.0098 (45%) 

Integration2 

1 

0.019 

0.007 (51%) 

0.0142 

0.0036 (64%) 

Integration2 

10% 0.006 

0.004 (21%) 

0.0028 

0.0013 (39%) 

Integrations 

1 

0.015 

0.006 (53%) 

0.0099 

0.0027 (66%) 

Integrations 

10% 0.005 

0.003 (23%) 

0.0018 

0.0007 (50%) 

Integrations 

1% 

0.003 

0.003 (<5%) 

0.0007 

0.0004 (20%) 


TABLE IV: Table of expected lower bounds on r with different 
integration experiments, following a stage of exploration. Re¬ 
sults are shown for three integration experiments with param¬ 
eters listed in Table [Ij assuming different levels of foreground 
residuals (from 100% to 1%) and of de-lensing {a = 1 for no- 
de-lensing or a = 0.1 for 90% de-lensing). The improvement 
achievable on average compared to the non-exploration case 
is shown in parentheses. 


cleanest patch with such an exploration, so that the up¬ 
per bound achievable on the tensor-to-scalar can be de¬ 
creased by as much as ^ 20% — 66% (or as much as 
a factor of 3 lower). We showed that adaptive-survey 
strategies for the exploration experiment can allow the 
cleanest patch to be identified a bit more efficiently if 
the experiment has the flexibility to rapidly change tar¬ 
gets. However, the improvement over a simpler rigid ex¬ 
ploration is not particularly significant. 

Although we have assumed in the analysis just one 
deep integration experiment, there are likely to be a sig¬ 
nificant number of independent projects, over the coming 
years, seeking inflationary B modes. The single explo¬ 
ration experiment to identify the cleanest patch is there¬ 
fore likely to benefit all subsequent experiments, allowing 
each of them ^ 20—66% improvements in sensitivity. The 
deep data from the exploratory experiment (particularly 
one at 353 GHz) should also provide valuable foreground 
templates for subsequent CMB experiments on any of 
the candidate patches explored. Finally, such data may 
also be useful to improve our theoretical understanding 
of dust-polarization by clarifying some of the open ques¬ 
tions and issues raised by Planck [7]. 

All of this improvement comes at the expense of just 
90 days of observation from the ground at 220 GHz or 
two weeks by a balloon operating at 353 GHz. This price 
seems affordable. 

We thank the anonymous referee for important com¬ 
ments which helped improve the clarity of our paper. 
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and NASA ATP grant NNX15AB18G. 
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